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USE OF SOURCE DISTRIBUTIONS FOR EVALUATING THEORETICAL AERODYNAMICS 
OF THIN FINITE WINGS AT SUPERSONIC SPEEDS 


By John C. Evvaed 


SUMMARY 

A series of ’publications on the source-distribution methods 
■for evaluating the aerodynamics of thin wings at supersonic 
speeds is summarised, extended, and unified. Included in 
the first part of this report are the derivations of: (a) the linear- 
ized partial-differential equation for unsteady flow at a sub- 
stantially constant Mach number; (b) the source-distribution 
solution for the perturbation-velocity potential that satisfies 
the boundary conditions of tangential flow at the surface and 
in the plane of the wing; and (c) the integral equation for 
determining the strength and the location of sources to describe 
the interaction effects (as represented by upwash ) of the bottom 
and top wing surfaces through the region between the finite 
wing boundary and the foremost Mach wave. 

The second part of the report deals with steady-state thin-nmng 
problems. The solution of the integral equation for upwash is 
presented and applied in numerous examples to obtain the 
velocity potential and loading for several families of curved 
plan boundary wings. The concept and the evaluation of the 
suction force along subsonic leading edges are included. The 
loading associated with the vortex sheet (or discontinuity in 
sidewash ) behind subsonic trailing edges is described and 
evaluated for wings satisfying the Kutta-Joukowski condition 
along subsonic trailing edges. 

The third part of the report approximates the integral equa- 
tion for unsteady upwash and includes a solution of the approxi- 
mate equation. Expressions are then derived to evaluate the 
load distributions for time-dependent finite-wing motions. 

INTRODUCTION 

The analysis of the aerodynamic effects in the vicinity of 
thin wings at supersonic speeds can be simplified by linear- 
izing the partial-differential equation of a compressible 
fluid. This linearization is accomplished by assuming that | 
the perturbation-velocity components associated with the 
wing are small in comparison to the velocity of sound and 
that the free-stream Mach number is sufficiently different 
from unity. Solutions to the linearized partial-differential 
equation of the flow that satisfy the boundary conditions 
in the vicinity of the wings must then be found. The flow 
must be tangent to the wing at its surface and no disturbances 
can propagate ahead of the forward Mach envelope of the 
wing. 

Numerous approaches lead to solutions of the steady-state 
thin-wing problems. For two-dimensional wings, the Ack- 


eret theory (reference 1) may be applied to evaluate both 
lift and wave-drag coefficients. The line-source and point- 
source methods of references 2 and 3 may be effectively used 
to obtain solutions for the aerodynamic coefficients of three- 
dimensional wings if the flow' on the bottom and top surfaces 
are independent. This condition holds for symmetric wings 
at zero angle of attack and for wing regions influenced only 
by supersonic leading edges. (That is, the component of the 
flow normal to the edge is supersonic.) Wing regions influ- 
enced by upwash in the flow field between the wing boundary 
and the foremost Mach lines, however, must be excluded. 

The effects of the upwash field may be evaluated for wings 
of straight-line plan boundaries by the conical-flow method 
of references 4 to 6. The doublet-distribution method of 
reference 7, as applied to the triangular wing, inherently 
includes the effect of upwash fields and may conceivably be 
extended to obtain quite general solutions for thin three- 
dimensional wings. An approach to this extension and a 
synthesis of the various methods as special cases of a general 
theory are included in references 8 and 9. The application 
of horseshoe vortices for solving the finite-wing problem is 
included in reference 10. 

The time-dependent aerodynamic effects (including the 
transient disturbances of gusts, changes in angle of attack, 
skin vibration, and flutter) near thin wings at supersonic 
speeds are more diffi cult to obtain. A number of investi- 
gators have studied two-dimensional time-dependent flows 
over thin wings; these flows are generally included as special 
cases of the theory of reference 11. The method of refer- 
ence 1 1 is similar to the steady-state method of reference 3 and 
includes three-dimensional or finite-wing solutions so long as 
the aerodynamic effects of the bottom and top wing surfaces 
are independent. 

If the upwash in the flow field between the wing boundary 
and the foremost Mach wave could be evaluated, the methods 
of references 3 and 11 could be applied to determine the 
aerodynamic coefficients in the flow field of lifting three- 
dimensional thin wings at supersonic speeds. This approach 
to the evaluation of the upwash field implicitly or explicitly has 
led to the methods of solution developed at the NACA Lewis 
laboratory and presented in references 12 to 19. The work 
included in these reports is extended and unified herein. 

The fundamental equations are derived in general form in 
the first part of the report. The applications of the integral 
equation for determining upwash to obtain the flow over 
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wing surfaces for time-independent and time-dependent 
problems are included in the second and third parts, respec- 
tively. Specific examples are included at pertinent points 
throughout the discussion. 

I— INTEGRAL EQUATION FOR UPWASH 

In order to unify the discussion, parts of the fundamental 
treatment of reference 11 are repeated. A similar discussion 
for steady-state potentials is given in reference 3. The 
presentation includes the derivations of: (1) the time- 
dependent, - linearized partial-differential equation for the 
perturbation-velocity potential of an ideal fluid, (2) the 
fundamental solution that will satisfy the boundary conditions 
on the wing, and (3) the integral equation for the upwash 
over the wing plan boundaries. 

Differential equation. — The linearized Euler’s equations 
for a compressible fluid may be written 
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where 

P 

t 

U 

x, V, * 

Po 

<p 


static pressure 
time 

free-stream velocity 

Cartesian coordinates (free stream parallel to x-axis) 
free-stream density 

perturbation-velocity potential (based on free-stream 
velocity) 

(For convenience, a complete list of symbols is included in 
appendix A.) 

When equations (1) are multiplied by dx, dy, and dz, 
respectively, added, and integrated, the result is 
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where g(t) is an integration constant at any given time. 
The linearized continuity equation assumes the form 
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equation (3) becomes 
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Elimination of ^ between equations (2) and (3a) gives 
Po 
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where M is the free-stream Mach number and the subscript 
has been dropped from c. 

Equation (4) is the required linearized partial-differential 
equation for steady or unsteady velocity potentials. If p and 
g are independent of ' time, the Prandtl-Glauert equation 
results. 

If the flow ahead of the body is uniform and undisturbed, 
the function g{t), which is independent of position, will be 

constant (from equation (2)) and equal to the term -fr 

Po 

will then contribute nothing to equation (4). Under these 
conditions, which are assumed for the rest of the analysis, 
equation (2) may be written 

f) « 

2 poU 

where C t is the pressure coefficient. 

A change of variable will convert equation (4) to a standard 
form of the wave equation. The transformations are 
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Equation (4) then becomes 
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Basic solutions of equation (6) corresponding to spherical 
waves are 


and 


where 




(7) 


R'=-Jx ,a +y'*+z f * 

(3 S =M ! -1 

/ arbitrary function representing a source strength 

Garrick and Rubinow have shown (reference 11) that the 
basic solution for M> 1 is obtained as the sum of equa- 
tions (7). If this solution is transformed to a general point 
in the x, y, z space, each of the addends becomes a retarded 
potential in the basic solution of equation (4), which is 
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where £, ij, and f are Cartesian coordinates of the source, and 
R is given by 

R=j(x-zy-!3\y- v y-p(z-l;y (9) 

Extended solutions of equation (8) may be obtained by 
integration with respect to any of the variables £, ij, or f. 
For thin wings, the sources of the disturbances will lie in the 
plane of the wing, which may be defined by the relation 
f=0. The immediate problem is to determine the strength 
of the source / in order to satisfy conditions of tangential 
flow on the surface of the wing, which in linearized theory 
is treated as though it were in the z = 0 plane. 

Fundamental solution for thin wings. — The impulses that 
the wing transfers to the fluid are primarily in the z-direction. 
This fact suggests that the strength of the source/ is closely 
associated with the z component of the perturbation velocity 

gj- The value of ^ at each local point on the constraining 

wing surface is determined only by the local wing slope near 
that point. 

The z component of the perturbation velocity obtained 
by differentiating the isolated source solution of equation (8) 
is 

b<p_ pzb<p . . 

^ (10) 

The wing dictates that ^ cannot be zero near z=0. Equa- 
tion (10), however, indicates that ^ will be zero in the 
x,y plane except at points near R= 0. (Infinite values of 
are generally included in this condition.) The quantity 

R is nearly zero either close to the source, where (r— £), 
(y—i j), and z are small, or on the boundary of the Mach 
cone from the source, where (at— %) i =P[(y— ^-fz 2 ]. As 
may be verified by direct calculation in specific examples, 

the second condition does not contribute to the value of ^ 

at the point (p,y, 0) on the constraining wing surface. 

In the vicinity of the point (a;,?/,0), the source strength 
per unit area assumes a mean value obtained by replacing 
£ and y by ar and y in the source strength / of equation (8) . 
The velocity potential at point (x,y,z) associated with this 
substantially constant source distribution in the z=0 plane 
may be obtained by integration over the source area included 
in the forward Mach cone. The field of integration (fig. 1) 
is bounded by the curves £= £i and (z— f) ! — p(y — ij)*— p *z * = 0 
yielding the velocity potential 

„ f MVS=*=»] dy 

* ' ' ,f( - X,y>t) h de /jyi V(*— 9*— /%—?)*— 0*2* 


(ii) 

For steady flows, this equation is the Ackeret value for the 
potential of an unswept two-dimensional flat plate and repre- 
sents the contribution of the local point to the velocity 


potential near that point. Partial differentiation of equa- 
tion (11) with respect to z gives 

^=w(x.,y ! t) = -2irf(x > y,t) (12) 

where w is the local z component of the perturbation velocity 
measured positively outward from the wing surface. The 
function/ at an arbitrary point (f, v) is thus proportional to 
the z component of the perturbation velocity at that point. 

The fact that ^ from equation (12) is independent of & illus- 
trates that the sources along the Mach cone for which R= 0 do 
not contribute to w. (Neither the basic solution of equation (8) 
nor its derivative, equation (10), is defined on the Mach cone. 
The integration, equation (11), of the basic source solution 
and the derivative, equation (12), of equation (11) are, how- 
ever, defined functions. The order of integration and differ- 
entiation of the basic source solution to obtain equation (12) 
cannot be interchanged without special precautions.) 

A condensation of the notation is advisable at this point. 
The time delays that occur in equation (8) can be denoted by 
r tt and rj where 

. fr-ea/ , ^{x-$y-f?{y-r,y-p? ' 

Ta pc ^ pc 

, y (i3) 

fr-QAf j(x-Q*-p(i/-y)*^p* 

T6 pc pc 

Functions of the time delays may then be abbreviated as 

fa=m,V,t-r a ) 1 

l (14) 

/c,6 =/ (£ l V, t — To) +/ (£ , V , t — Ti) J 

In this notation, the velocity potential that is obtained by 
integrating the sources in the z=0 plane over the area S 
included in the forward Mach cone is 

= 1_ f P (w a -\-Wt)d^ dy 

2t J Js Vfc ~ Q l ~P (y—y) l —pz 2 

_ l r r Wa -» dj dr > rtK\ 

2 ^J Js V(*-0 1 ; 

Equation (15) was derived in reference 11 for time-dependent 
velocity potentials and in reference 3 for time-independent 
velocity potentials. 


Forward Mach cane 


* or § 

y 

U z=Q plans (e quo lion R=0) 

, or £ ' 

Fig use X. — Ffeld of integration for evaluating velocity potential (equation (10)) near plane 

of sources. 

A physical interpretation (reference 11) may be given for 
the time delays, r a and t„ of equations (8), (13), and (15). 
If a disturbance is generated at point (|,jf) only at time f=0. 
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the wave front from that disturbance will be spherical about 
a center that moves with the free-stream velocity. The 
trace of the waves on the 2=0 plane is illustrated in figure 2. 
The wave front will enter and emerge from the point (x, y, z ) 
at two later times r a and r b . The equation for the spherical 
wave path that passes through the point (x, y, z) is 

(x— £— L7t)*+ (y— ij) a +s 2 =cV 

or 

= {x— £)M V (a— -g) 2 — ft 2 (y— q) 2 — ffV 
r pc = pc 


tribute to the velocity potential. These two waves originated 
at (£,?j,0) at times (t— t„) and (f— tj). 

The potential on the surface of the wing may be evaluated 
by equation (15) with z set equal to zero. The boundaries of 
integration then include the forward Mach lines from the 
point (*,y,0) . The integrations may generally be simplified 
in an oblique set of coordinates whose axes lie parallel to 
these. Mach lines. The transformation equations relating 
the oblique and Cartesian coordinate systems shown in 
figure 3 are 


which agrees with equation (13). Thus, at a given point 
(x, y, z), the strength of the same wave at the two times r a 
and (in accordance with equation (15)) contributes equally 
to the velocity potential despite the change in the radius of 
the wave front. At a given time t, only the wave fronts 
that are entering and emerging from the point (x, y, z) con- 



Fiqcre 2 .— Relation between time delays t. and t,, and position of wave front. 



Figube 8.— Comparison of Cartesian and oblique coordinate systems. 
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(In references 12 to 19, the symbols u and v were employed to 
represent the oblique coordinates. These symbols were re- 
placed by r and s in the present report to avoid confusion 
with the common usage of u and v as velocity components.) 

25 

The elemental area in the r, s coordinate system is jp dr ds. 
Equations (15) and (13) then become (for 2=0) 


1 C C w a-t dr ds 
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If only supersonic leading edges are included in the 
forward Mach cone from point {r w ,s w ), w may be evaluated 
in terms of the effective wing slopes <r measured in ij= con- 
stant planes by the relation 


w=U<r 


(19) 


(The value of a would thus be positive for either top or 
bottom surface of a wedge wing at zero angle of attack.) If 
a subsonic leading (or trailing) edge is also included in the 

V} 

forward Mach cone, the slopes of the streamlines 

associated with the upwash between the wing boundary and 
the foremost Mach wave must be evaluated and included in 
the calculation (equation (17)) for the velocity potential. 
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Integral equation for upwash between wing boundary 
and foremost Mach. line. — Equation (17) contains no 
description of the origin of w, the velocity potential in the 
flow field S D of figure 4 may then be independently calculated 
with respect, to either the top or bottom wing surfaces. 
These two potentials are 


U f C (<rr)g,& dr ds 
<pT 21frrJ Js. V (r D —r) (s D —s) 

U r r \g, » dr ds 
2M~ J Jsoi/(r D — r) (s D —s) 

_ V C C (ctb) a, b dr ds 

' fB 2 A/xJ Js. V( r D —r) (s D —s) 

U C C K,t dr ds 

2Mir J Jsd Vfa,— r) ( s D —s ) 


( 20 ) 


( 21 ) 


where r D and Sd are the coordinates of the point at which <? 
is evaluated, <r T and a s are the slopes of the wing on its top 
and bottom surfaces, and \ represents the slopes of the 
streamlines in the field Sd from the point of view of the top 
wing surface. The pressure may be calculated by sub- 
stituting' either equation (20) or (21) into equation (2). 
Because no pressure discontinuity can persist- across the S D 
region of the z= 0 plane, 

r I T t dfffl , Tr da s 

dt + u dx~ dt ±u dx 
Equation (22) has the solution 

<Pt=<Pb+2H{x— Ut,y ) 

where H is an integration function. From equations (20), 
(21), and (22a), there result 
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The function 2 H represents the difference in potential 
across the 2=0 plane (equation (22a)) corresponding to the 
strength of vorticity in the wake of the wing. For anti- 
symmetric wings, <r B -hcrr=0 and H is the potential on the 
top surface of the vortex sheet. 

The foremost Mach wave (fig. 5) originating on the leading 
edge generally represents a line of infinitesimal disturbance 


(20a) 

(21a) 

(22b) 


( 22 ) 

(22a) 


along which H may be set equal to zero at all times. The 
function H remains zero along y = constant lines for values 
of x not intercepted by the wing or material body (region 
Sd , 2 of fig. 5 (a)). The region Sd , i generally contains a 
vortex sheet lying in the plane of the wing, and II is not 
zero. The function H{y), established along the wing edge 
at some time t, remains unaltered for later times along a 
curve that sweeps downstream with the free-stream velocity 
and has the form of the wing trailing edge. 



(a) Fields of integration for evaluating upwash between wing boundary and foremost Mach 

wave. 



(b) Fields of integration for equations (28) and (29). 

Figuee 4. — Wing region influenced by isolated subsonic leading edge. 


Equation (22b) represents the integral equation for evalu- 
ating the slopes of the streamlines in the upwash field. These 
slopes may generally be defined by the following two equa- 
tions: 



(Xp) a , a dr ds 
■y!{r D —r){s D —s) 
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The actual slope X is then the sum of the solutions for the 
slopes obtained from equations (23) and (24) : 


X — X 0 +Xa . 


(25) 


For regions such as S D of figure 4 (a), 27=0 and only equation 
(23) need be considered. For steady-state solutions, the time- 
delay notation a, b may be disregarded in equation (23). 

II— TIME-INDEPENDENT FLOWS . 

The concept leading to the derivation of equations (23) 
and (24) may be extended in principle to define the upwash 
velocity components in the plane of arbitrary thin wings at 
supersonic speeds. Equation (17) then yields the velocity 
potential from which the aerodynamic coefficients may be 
derived. The effects of the upwash field have been evalu- 
ated for a class of curved plan-form wings, and the corre- 
sponding velocity potentials of time-independent configura- 
tions have been derived. Included in the calculations are 
wings with subsonic leading and trailing edges connected to 
supersonic leading edges and a few cases of wing regions 
influenced by interacting upwash fields. 

Evaluation of isolated regions of upwash off subsonic 
leading edges. — The steady-state solution of equation (23) 
is conveniently derived for the wing of figure 4. The origin 
of coordinates is placed at the junction of the supersonic and 
subsonic leading edges. These are defined, respectively, by 
the equations 


. s=Si(r) or 

s=s 2 (r) or 

Equation (23) then becomes 


r=ri(s)] 
r=r 2 (s) ! 


(26) 


f r ° dr X ds po dr f 
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«i(r> 2 tJs d —s 
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(Inasmuch as is zero, X 0 is replaced by X in equation (27) .) 

Because equation (27) must hold for all values of r D , which 
does not appear in the integration with respect to s, the 
equation may be reduced to 


C’ D X ds C 
J*!(r) -yJs D — S J»: 


X ds (tj B — (7j ,)ds 

»itr) 2^S d — S 


(27a) 


This reduction shows that for steady-state solutions only the 
wing slopes along the forward Mach line contribute to the 
upwash at the point (r D) s B ). 

Equation (27a) may be solved, as in reference 14, as an 
example of Abel’s equation (references 20 and 21) to give 


X(r,s) 


___1 ^ ds D C * iW (?b — g>) do 


k ds J W) ■Js—SpjW) 2 Vsfl— a 


(27b) 





to) 

(c) Point (r„, «„> on wing surface. 

Figure 5.— Integration boundaries for evaluating wing region Influenced by subsonla trailing 

edge. 
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where v has replaced s as a variable of integration. Equa- 
tion (27b) may also be expressed in differential form as 


d\ (<Tg — & t) V 5a — ® 

do 2ir (s—v)tJs —82 


(27c) 


where a wing element of length do and having the effective 
slope < - B 2 Ct at point (r,v) produces an increment d\ at point 

(r,s) (fig. 6). Either by interchange of the order of integra- 
tion in equation (27b) or by integration of equation (27c), 

an alternative expression for X results 

\(r,s) = f *’ i/r»—v T )^at—vdo (27d) 

2irVs — s 2 J*i s — V 



Floral 6.— Geometric significance of actors in equation (27d). 


If (o* b ' tr x) is independent of s, either equation (27b) or (27d) 
gives 


\(y) - -tan -'Jh® W.1 (27e) 

iT L V #— *i(r) V s— * s (r) . 


Aerodynamic effects of isolated regions of upwash. — 
The contribution <p D of the upwash field S D to the velocity 
potential on the surface of the wing in figure 4(b) may be 
computed from equations (17) and (19) as 

_ U P r 2CO dr f*» \ds , . 

‘ pD ^r a —rjH<f) -Js a — 8 


If X is eliminated from equation (28) by means of either 
equation (27a) or equation (27d), there results 


U_ dr p>2« (g- g — g- T ) ds 

‘ pD Mvj t -yjr w — rJ*iC) 2 -yjs u —s 

U C C (o~g— <r T ) dr ds 

MtJ Js oJ 2 V 0»— r) (Sm — s) 
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The contribution of the upwash field is thus replaced by an 
equivalent contribution on the wing surface. The surface 
velocity potential at (r„, s w ) then becomes 

_ U f C or dr ds 

9 Js„. 1 V(r B — r)(s„— ») 

U C C (q-g+q-r) drds 
SfirJ J s ®j 2 ^(r a — r ) (s M — a) 


.For thin-plate wings, or for antisymmetric solutions (such 
as evaluation of the effect of unif orm rate of roll), (q-g+qr)=0 
so that the integration over the area S„.a vanishes. In this 
case, the contribution to <p of the upwash field cancels the 
direct contribution of the part of the wing that generated 
the upwash field. Hence, alteration of the thin-plate-wing 
plan boundaries in the region S a-2 (fig. 4 (b)) does not alter 
the surface velocity potential at point (r „,.?„) . (An elegant 
derivation of this fact is presented in reference 22. Because 
there is no contribution to the velocity potential ahead of 
the forward Mach cone, equation (17) may be written for 
steady boundary conditions as 

1_ C r ‘ p- w dr ds 

MtJ -c J _<= y (r„— r) ( s u —s ) 

Two auxiliary point functions may now be defined as 


X(r a , s) = — 


F(r, «*) = — 


1 1 

(* r ‘ w dr 

Mir] 

Vr.-r 

1 1 

[**■ w ds 

375 rj 

— \'s a —s 


In terms of these point functions, the velocity potential is 



Xds f r - 

7= or 

y S jo 8 J — 00 


7 dr 
v / JV zr r 


Each of these equations for X, F, and <a are forms of Abel’s 
equation and may be solved accordingly. For an anti- 
symmetric wing of plan boundary as shown on fig. 4 (b), 
<P (r, s„)=0 if r<rj(8„). Hence Y is zero for this range also. 
The velocity potential may then be evaluated only from the 
region of fig. 4 (b).) Pressure coefficients may be ob- 
tained by substituting <p from equation (29) into equation 
(2a). Perturbation-velocity components are obtained from 
the gradient of <?. 

As illustrations of equation (29), two examples are pre- 
sented. The first is a wedge wing (fig. 7) with plan-boundary 
equations s——k x r arid s=k 2 r. The wing slopes are taken 
as q-g=q+a and a T =cr—a, where 2<r is the wedge angle and 
a is the angle of attack. The velocity potential on the top 
surface of the wing is then (reference 12) 


2 U<r\'s a -\-k l r u ^ lEfl , (s a — hrj) w 

—tg [_N*r t! " ‘Vtt +— x 
lo& VE+ vi] + [ V T? J + 

(Sg-ykird) [ki{k 2 r a — 8ar)1 

vfei V (*i+* 2 )s»J 


(30) 


This solution for <P is the addition of terms that are dependent 
on and independent of the angle of attack. The effects of angle 
of attack may therefore be evaluated on flat-plate wings 
whose plan boundaries coincide with those of the chosen 
wing. This example illustrates a general superposition 
principle that follows directly from the linearity of equa- 
tion (15) with respect to either w or a (from equation (19)). 
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The velocity potential of a flat-plate wing (fig. 8} having a 
straight-line leading edge s=— hr coupled to an arbitrary 
subsonic leading edge s=s 2 (r) or r=r a (s ) may also be 
evaluated explicitly. In this example, £r s = — <rr=a, so that 
</> is 

[Vft.+l.rO *V-«J 

where r % is a function of s a . The factor multiplying a in 
equation (30) may also be obtained from equation (31) by 
s 

setting r 2 —jr‘ 

Some complications occur ilthe supersonic leading edge is 
connected to two subsonic leading edges that are located in 
the forward Mach cone from the point (r w , s w ) (fig. 9). The 
contributions to the velocity potential of each of the upwash 
fields may be independently calculated from equations (28) 
and (28a) to give - ■ 

U_ P P X i dr ds U P P (o~ fl — g r ) dr ds 

Mir] Jsd.i V (r„ — r) n«> 2-\l(r a —r)(s u> —s) 

(32) 

U P P Xi dr ds _ U P P Q <r B —<r T )drds 

Mvj Js 0 ,t ^{r a —r){s w —s) MirJ 2 V(r»— r)(s»— «) 

-"i 

■ ■ (33) 



V or p 



Section A -A 


Figure 7.— Dlscontinuously swept wedge wing. 



I" 



Figure 9.— Regions of Integration for calculating velocity potential on surfaco of wing of 
finite thickness Influenced by two independent perturbed flow fields external to wing 
surface. 


Addition of these contributions to those of the wing area 
SaKi+s+s+D, as evaluated by equations (17) and (19), gives 
(reference 12) 


U P P (g'g-l~o'r)dr ds 
’MirJ Js. (i+i) 2 V (r«— r) (s w — s) 

U P P <sb dr ds 

^ J Js,.> V(r„ — r) (s„ — s) 

U P P <tt dr ds 
Mir j ij(r a —r) (s^—s) 


(34) 


The role of <r T and <r B is thus effectively interchanged for the 
regions of integration S a ,t. Likewise, for antisymmetric con- 
figurations (o-B+irr=0), the areas S Wt i and S w>i give no 
contribution to the velocity potential. 
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Aerodynamic effects of interacting regions of upwash. — 
If interacting upwash fields are included in the forward 
Mach cone from the point (r u , s a ), as in the shaded region of 
figure 9, parts of the upwash field must be evaluated before 
the velocity potential can be computed. The effective wing 
in application of equation (34) includes the upwash fields 
S D 'i and S D 2 (fig. 10) for which ub=—\. The velocity 
potential may then be written * 


U C r o T drds , 

Mir J Js„. 4 V( r »— r)(8w— 8J 

U C C X t dr ds . 

Mir J Jsdj VOv— r) (s„— s) ' 

U f f* X 2 dr ds 

J Jsoj -J(r u — r) ( s a —s ) 


(35) 


(o-fl-f-gr) drds 


= u r r 

1 Mir J Js.o+n 2 V ( r ii— r ) (s m —s) 

U r C <t b dr ds 

Mir J Js.j V( r w —r){s u —s) 


Each of the upwash streamline slopes may be evaluated by 
equations similar to equation (27d). In the notation of fig- 
ure 10 (and with e replacing s as a variable of integration), 
the contribution of the flow field S D ,i is 


U 


Mir., 


f 

\i dr do 

U fr 2 (n) dr fi 

JSd 

r)(s u — 0 ) 

Mir Jo ^r u -rjt-. 


* 


U f n(n) dr f*t( rj do (ff a — (Tr) V s a — s ds 

2MV 2 Jo ^{r u — r) V ( s » — °) (P~ 82) J*i(d (p—a) 


(36) 


Interchange of the order of integration with respect to c and s and evaluation of the integral with respect to o give 

U_ f f X t dr ds U_ f f ( S s~ g ) ( S *+S 2 — 2s) — 2(s 2 — a) (s„— s) ~j (<x B —<r T )dr ds nR ^ 

Mrj Jso.t Tl(r a —r)(s u —s) 2Mir 2 J Js„.,iL 2 ' 2 V(« 2 — s) (s m —s) (s*— s 3 ) (s 3 — s 2 ) J V (r a —r) (s a —s) 


where s 3 is a function of r u and s 2 is a function of r. (Equation (36a) reduces to the form of equation (28) when s w =s 3 (rj).) 
Likewise, 

U_ C f X 2 dr ds U_ f f fir , x .__, (r 2 -r) (r„+r z — 2r) —2 (r„-r) (r a -r) 1 (<r B -<r T )drds 

MirJ J Soj V^io— ?■)(«*— s) 2Mir J J Js..nL2’ r 2 V(r s ~r) (r„-r) (r B -r^ (r 2 -r s ) J V (r*—r) (s B —s) ‘ 7 


where r 2 is a function of s„ and r 3 is a function of a. Elimination of Xi and X 2 between equations (35), (36a), and (37) yields 
the velocity potential in terms of the wing slopes and plan-boundary equations. An alternative scheme (suggested in refer- 
ence 14) would be to derive all of the upwash slopes X by means of equations similar to equations (27b) or (27d) and then to 
evaluate a directly from equations (17) and (19). More complicated examples including doubly interacting external 
upwash fields may be solved in a similar manner. 


If the upwash fields continuously interact as they do for 
the wiog shown in figure 11, the solution of an integral 
equation is required (reference 14). The slopes of the 
streamlines from equation (27d) in the upwash region 
Sd.i are 


Xi(r,s) 


1 r 

IT Vs — 82 


(da — <r r) V®2 0 j" 

2 ( 8 - 0 ) a 



X 2 (r,o) Vs 2 — 0 
(8—0) 



( 38 ) 


where 0 is the integration variable in the s direction and s 2 
and s 3 are functions of r. 

Similarly, 


X 2 (r,s) 


-1 rr 

irVr— r% l_J ri 


(p B — o- r ) Vr 3 — u 
2 {r—u) 


du— 



\(u,s) Vr 3 — u 
{r—u) 



(39) 


where u is the integration variable in the r direction and r 2 
and r 3 are functions of s. The functions X r and X 2 can there- 
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fore be determined from the two equations, at least in 
principle, either by successive approximations or by direct 
substitution and solution of the resulting integral equation. 
(If the wing has a symmetrical plan form about the a-axis, 
X 2 (r,s)— ^(s, r) and the two equations (35) and (36) become 
the same. If the flow is conical, either \i or X 2 is a function 
of s/r.) Once the slopes of the streamlines in the upwash 
field are determined, the velocity potential may be evaluated 
from equations (17) and (19). 

Calculation of perturbation-velocity components. — The 
perturbation-velocity components may be obtained by 
partial differentiation of the velocity potential. This differ- 
entiation may be performed prior to the integrations indi- 
cated in equations such as (17), (29), (34), and so forth. 
The method is illustrated, as in reference 15, for the wing 
of figure 12, which includes only supersonic leading edges in 
the forward Mach cone. If the velocity potential is computed 
at point (x,y) and at point (z-f dx,y), the result is 



where the line integral is evaluated along the line abd, and 
the coordinate £ may be eliminated from the line integral 
by means of the equation for the leading edge. Subtraction 
of equation (40) from equation (41) and cancellation of dx 
yield 

b<p_ u r r 5 ?^ 

dx * j Js V(z-Q s -/3 2 (y-7) 2 


U r a dy 

irjabd YCk— £) a — j 8 2 (y— 17) 2 


(42) 


This result may also be obtained by formal differentiation of 
equation (40). In a similar manner 


U f C , U f <r# 

Sy ^ J J s V Jabd Y(z— £) 2 — jS 2 (3/ — 77) s 

(43) 


Equations (42) and (43) may also be written in the oblique 
coordinates as 



(43a) 


I* 



Fioure 11.— Wing plan form and areas of upwash for equations (38) and (30). 



M 

(a) Wing point Or, »). 



x+dx,y 

(b) 


fb} Wing point (z+tOr, y). 

Figure 12.— Fields of integration for equation* (40) to (42). 
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The line integral of equation (42) corresponds to the discontinuous change in flow direction that occurs along the 
leading edge, and the surface integral corresponds to continuous changes in slopes over the wing plan area. (Evaluation 

of the surface integral along lines of discontinuous ^ results in additional line integrals.) The impulses that lead to the 

pressure distributions on supersonic wings are thus generated at points where the flow direction is altered. For flat surface 
wings, these impulses are generated only along the leading edge. 

When the surface velocity potential is influenced by an upwash field, as in figure 4, the velocity potential may be 
calculated from equation (29). The evaluation of the perturbation-velocity components is similar to the derivations of 
equations (42) and (43) except that there is effectively a shifting line of discontinuous wing slopes between the plan areas 
S W 'i and Su'i, which adds line integrals to the calculation (reference 15). The partial differentiations of <p with respect to 
x and y (fig. 4 (b)) have been evaluated in appendix B and are expressed in oblique coordinates as 


U rr (£+&)** 
<dr 2/3rJJs„.i V(r„— r) (s„— a) 


U 

4/J 


P (<r fl +<r r ) , &(e- a +<*>)! j- 

S r + ds \ ar(tS U C a T (ds—dr) 
vJJs.., V (r„— r) (s u — s) 20irJ„b V(r»— r) (s a —s) 


U_ 

4/3 !T, 


r (<r B -j-OT 
JbOd V ( r »— 


-f— g 2 >) (i ds—dr ) U 


(l-—') 

V dig/ C (a T —ff B )ds 
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dy 2irJ Js.., 

£f 

2irjb 


r) (s a ~s) 4/Sir 

(£-£)»* 


VSu — 8 


(<r a -f-o~r) 
bod 2 


7 dsj r 
Vr K — r 2 Jb 

dfas + cr) SfcTfl-j-fTr) - ! j. 

r r L £>r J ar as n - ffr ( ds+dr ) 

J Js*,i Y(r„— r) ($u — 8 ) 2?rJ a b Y(r tt — r) (s ffi - 

< d I±) 

‘ dSg/ f (cr r — cr B )ds 

V(r»— r) (s„— s) 1_ 4 :t ^jr a — r 2 Jbd V s »— * 


(44) 


V (r„-r) (s„— s) 4x w 


(rfs+dr) 


(45) 


where r s is the wing-plan-boundary equation r=r 2 (s„). The 
more complicated expressions, such as equation (34) and 
combinations of equations (35), (36), and (37), can be 
treated in a similar manner. 

The expressions for the perturbation-velocity components 
are simplified for flat-plate wings having slopes o- B = — <r r =a. 
(These components may also be regarded as the contribution 
associated with angle of attack on a cambered wing.) Equa- 
tions (44) and (45) may then be written 


5oj_ 

Ua ( 

’ ( ds—dr ) 

Uaf 

£ 

dr 

20ir J 

ab V : r) (S»— ») " + 

’fa V 

1 

and 
d < p 

Ua 

C ( ds+dr ) 

Ua, 

^ ds „ 4 

by~ 

2ir, 

Jab V(r„— r) (a„— s) 



g — 8i(r a ) 

r u — r 2 


(46) 


(47) 

The line integrals of equations (46) and (47) may be numeri- 
cally or graphically evaluated by the procedures outlined in 
reference 16. If the supersonic leading edge conforms to the 
equation s——kir, equations (46) and (47) give 


^p_Ua r /. dr 2 \ Is 

~ p* Lv 


dr 


*i Tj 

ds w J \ r w—r 2 


(*i+l) 


VS 


tan" 


j ki(r u -r^ ~ 

V s»+*ir 2 _ 


a v5__. U<*Vfi A-dri\ hu+kirt Ah~l) /*i 0 

&iT *■ LV + dsj \i 7^7T + 


(48) 

+hr 2 J 

(49) 

Equations (48) and (49) could also be obtained by partial 
differentiation of equation (31). 

Along the subsonic leading edge, (r„— r 2 )=0 and an infinity 
occurs in the perturbation-velocity component of equations 
(46) and (47). This infinity represents a breakdown in the 
validity of the linearized theory, and corresponds to the for- 


mation of a stagnation point near the leading edge of two 
dimensional subsonic airfoils. Even though unbounded pres 
sures may be obtained from equations (2a) and (46), these 
values should be interpreted as representing some fraction of 
free-stream stagnation pressure (or near-vacuum conditions). 
As was observed in reference 23, discontinuities in wing-plan- 
boundary tangents along subsonic leading edges lead to dis- 
continuities in the x and y pertubation-velocity components 
and in the wing loading along the Mach line from the dis- 
continuity. 

A physical consideration of pressure disturbances from the 
wing-plan boundaries indicates that the velocity component 
lying parallel to the edge must be continuous across the wing 
edge. This concept is evident from oblique-shock relations 
for supersonic leading and trailing edges and may be verified 
by equations (20), (46), and (47) for a class of subsonic lead- 
ing edges. All changes in the velocity at the wing-plan 
boundary must therefore occur in the velocity component 
perpendicular to the edge. 

Evaluation of suction forces on subsonic leading edges. — 
Near the wing edge, the line integrals of equations (46) and 
(47) are zero. Likewise, the direction slope of the perturba- 
tion velocity is 


dj/ = P \ 1 dsj _ _ A dx\ 

(i dri \ wA 

dr v dSu) 


(50) 


where 


fdx\ 

WA 


is evaluated from the wing-plan-boundary 

equation. The perturbation velocity is thus normal to the 
wing-tip boundary and has the value 

[0 -&M 1 +U] (61 > 

r »-* r s 
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In accordance with references 24, 7, "and 19, this type of 
variation in as n-> 0 leads to a suction force per unit 
length along the leading edge of magnitude 

dF n _^pU 2 o?(s w —s l ) 


Urn /gr, 

dl ttM \ ds a 


(52) 


where is the force normal to the leading edge and dl is an 
infinitesimal distance along that edge. If dl is expressed in 
terms of ds and F is defined as the suction force in the flight 
direction, equation (52) may alternatively be written 


lira 

s w ~* 9 2 


dF 


3s' 


4:pU 2 CX 1 (s w — S]) 
ttAF 


(’-ISX/I (53) 


The suction force along subsonic leading edges can thus be 
evaluated explicitly by equation (53) for wing boundaries 
that are not influenced by interacting upwash flow fields. 
Equation (53) has been applied in reference 19 to demon- 
strate that appropriately curved plan boundaries may lead 
to wings of higher lift-drag ratio than do straight-line plan 
boundaries. 

Velocity potentials associated with discontinuous side- 
wash. — If the upwash field lies in the wake of the wing (fig. 
5 (a)), the integration function II (y) of equations (20a), 
(21a), and (22a) will not generally be zero (reference 15). 
A discontinuity in sidewash will then persist across the 
z— 0 plane corresponding to the strength of the vortex sheet 
trailing behind the wing. The function II may be adjusted, 
if desired, to obtain solutions for the velocity potential that 
will satisfy the Kutta-Joukowski condition along subsonic 
trailing edges. This condition requires that for flat-plate 
wings the perturbation-velocity components are continuous 
across the trailing-edge boundary. For finite-thickness 
wings the velocity components must be finite at the subsonic 
trailing edge. The immediate problem is to. establish the 
part of the velocity potential on the surface of the wing that 
is associated with the function H(y). 

The region S D .i of figure 5 that contains the vortex sheet 
may be temporarily considered as a portion of the wing. 
The virtual wing tip is then the y= constant line denoting 
the junction between the leading and trailing edges. If 
equation (29) is applied to the virtual wing-tip region S Di i 
of figure 5 (a), there results (in the notation of fig. 5 (b)) 


_ _ u r r 
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(o ~ b ~j- v r) dr ds 
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T C C (Xp+X g )dr ds 
Mr J Js D4 (s D —s) 


(54) 


• If this expression is compared with equation (20a), the 
result is 

J7 ^_J7T C f (Xp-f X g )dr ds C C (<T B —<r T )dr ds ~1 

Mr [_J Js Cil -)\r D —r)(s D —s) J JVi 2 V(r c —r)(s fl — s) J 


( 55 ) 


In the notation of figure 5 (b), equation (55) may be written 
'*n> dr r f**(b (o> — c B )ds 
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+ 


(56) 


By use of Abel’s solution (references 14, 20, or 21), equa- 
tion (56) may be inverted to yield 

(vr — r B )ds . r*o (Xo+X g )d8 ~| 
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(56a) 


where H(y) has been expressed in oblique coordinates. 
Equation (56a) may be integrated by parts to give an alterna- 
tive form 
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(56b) 


The velocity potential at points on the wing influenced by 
the subsonic trailing edge (fig. 5(c)) is 
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-(57) 

The second member of equation (56a) or (56b) may replace 
the first member along fines of constant s D (or s B ) that 
extend, across the wing. Equation (57) then becomes 

= i rvo dr fb p n (^JT^ dr ° \- 
9 v J'u- dr a —r LdrJ»„ s!r—r D J 

U p(Q dr p 
Mr Jr -yjr a —r )h 


U p» dr p» a T ds 
Mr -yjr w — rj*i(0 -yjs a — 


hi 7 ) (q-g-j-g-r) ds 
W 2 y/Su—8 

ds 


(58) 


This expression reduces to equation (29) if Id is zero. Inas- 
much as II may be arbitrarily chosen, an infinity of solutions 
for regions influenced by subsonic trailing edges can satisfy 
the boundary conditions for thin wings at supersonic speeds. 
On the other hand, if the Kutta-Joukowski condition is ap- 
plied, only one solution remains — the velocity potential that 
allows the flow to leave the subsonic trailing edge smoothly. 
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Once tlie function H is chosen, equation (56a) may be 
applied to evaluate the upwash (or X) over the wing edge. 
The slope X has been considered as the sum of two sets of 
slopes, one independent and one associated with H. 

X=Xo _ bXa' (25) 


Equation (56b) may then be written as separate equations 
for Xo and X H as 


and 
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(27a) 


(59) 


The bracket of equation (57) could thus be replaced by the 
first member of equation (59) if desired.* 

The solution of equation (27a) is obtained from equa- 
tion (27d). The solution of equation (59) follows from 
Abel’s relation and may be written in the form 


X*= — 


1 [ b 
Uir bs 



(60) 


This is the steady-state solution of equation (24), which could 
have been obtained directly. If H(0) =0, equation (60) 
may also be written as 


bH , 
Nir dr D 


X*=-#£ p (60a) 

Lit bs Jj 2 (r) »£> ^r—r D 


Substitution of X 0 from equation (27d) and \ s from equa- 
tion (60) into equation (25) yields the slopes of the upwash 
streamlines in the wake of the wing (regions'll of fig. 5 (a)). 
If the region S D ,i is treated as part of the wing, the slopes of 
the streamlines in the region S Dy2 may be evaluated from 
equation (27d). As in the derivation of equations (36a) 
and (37), the influence of interacting upwash fields may 
then be determined. 

The velocity potential associated with H in equation (58) 
may be written in either of the alternative forms 


1 (*,(.-) S |> Hdr, 

r dr Vr— 
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where s a has replaced s D as in equation (58) . The question 
arises as to the nature of the velocity potential and the 
perturbation-velocity components that result from the 
function H, which depends only on the coordinate y. Two 
approaches are followed to answer this question. The first 
approach is to evaluate Vg when H assumes the form of a 
power series. (This approach is useful for explicit evaluation 
of <p for wings with polynomial plan-boundary equations.) 
The second approach is accomplished by direct integration 
and differentiation without further assumptions on the form 
of H. 

The function H (y) =H ) is assumed expansible 

in a Taylor’s series of the form 


H 


( 


S u T j) \ , Cti(Sa, Tg) , @2 F d) 2 i 

M / a ° + M ‘ A/ 2 


(62) 


(Because H(0) is zero, a a will be zero for isolated wings.) 
Substitution of H from equation (62) into equation (61) 
yields for the velocity potential (see reference 15 for details) 



(63) 


where F represents the gamma function, o is an integration 
variable, and r 2 is a function of s w . The integrations of equa- 
tion (62) may be expressed as incomplete beta functions if 
desired. 

The x component of the perturbation velocity that is 
associated with H may be obtained by differentiation of 
equation (63) as follows: 



The infin ite, series is a function only of s w for a given plan 
boundary and Mach number. The coefficients a„ control 
the strength of the vorticity in the wake of the subsonic 



trailing edge. The factor — w =- represents an inverse 

Vr»— r 2 

square-root singularity along the wing-plan boundary that 
leads to an infinite pressure coefficient in the linearized 
theory. 

The x component of the perturbation velocity in regions 
influenced by a subsonic trailing edge (fig. 5 (c)) is obtained 
as the sum of equations (64) and (44) (with ec 0 set equal to 
zero) 
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For flat-plate wings, <J B — — O' T — CC and equation (65) becomes 
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(66a) 


JJct 

where g(s„) represents times the bracketed part of equa- 


tion (66). The quantity ^ 1 is constant along lines 

of constant s„, so that the strength of the shed vorticity 
could feasibly be experimentally determined by a single 
pressure measurement on the wing surface for each value 
of s w . 


If the Kutta-Joukowski condition is to be satisfied for 
antisymmetric wings, all three of the perturbation-velocity 
components must bo continuous across the subsonic trailing 
edge. Continuity in any one of the three components will 
evaluate the coefficients a„ and assure the continuity of the 

other two. Continuity of ^ across the subsonic trailing 

edge for flat-plate wings, covered by equation (66a) for 
example, requires that g(s w )=0. The solution that satis- 
fies the Kutta-Joukowski condition is thus 
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The coefficients a u a 2 , . . . a n may be evaluated by the 
relation g(s w )= 0, or for flat-plate wings, 
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r (n+~) M n 


( 68 ) 


where p— ^/r 2 —s a . Equation (68) is a power series in p 
whose nth coefficient a n may be evaluated in terms of tho 
(2n— l)th derivative about p= 0. Because the first member 
of equation (68) is finite when p—0, a c —0 as previously 
stated. Differentiation of the first member of equation (68) 
can be accomplished by successive applications of the relation 


d _ 2j t —s w 

dp fdr 2 _ \ 
\ds a ) 


d 

ds w 


(69) 


By application of equation (69) and L’Hospital’s rule to 
equation (68), the coefficients a 2 and Oj, for example, are 


given as 

j l _ds i dr t 
/ dr 2 ~3J 

/ dr t 


s’ 

1 

II 

iH 

13 

(70) 

V 

evaluated at’r 2 =s ti =0. 

— 1 

(LSy} 



2 \ dvi ds w J dsj L\uSa> / \«s K / 8 \ds a /J 

MU a ( dsi drt\-$ pggi (d£j\, dgi d t r i ~\ V( dr t ( dr 2 3 ds^ dr 2 \~ | , . 

2 \ dr 2 dsj L^ r s 2 \dsj ~^dr 2 ds^J [_\^ L ) \ds w / 8 \ Hr 2 dsj J ' ’ 

evaluated at r 2 =s w =0. The other coefficients may be similarly evaluated. If wings with straight-line plan boundaries are 
considered (as for conical flow wings), all the coefficients a u a 2 , a 3) . . . a n except the first are zero. For curved plan bound- 
aries, higher-order terms in the series expansion for H are generally required to satisfy the Kutta-Joukowski condition. Tho 
expressions (70) and (71), which determine the first two coefficients in the series expansion (equation (62)), are unbounded if 

approaches unity at the origin. 

Once the coefficients a t , a 2> a S) . ,. . a n have been determined, the velocity potential and its partial derivative with 
respect to x are given by equations (63) and (64), respectively. The y component of the perturbation velocity is obtained 
by partial differentiation of equation (63). 
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The 2 component of the perturbation velocity associated with the function H may be obtained by substituting equa- 
tion (62) in equation (60) (a 0 =O): 
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(72) 


The terms factored by a a , g 3 , g 4 , . . . g„ do not appear for 
wings with straight-line plan boundaries. 

An alternative approach to the determination of the form 
for the perturbation-velocity components associated with the 
function II may be obtained by integration and differentia- 
tion of equation (61a). This manipulation is accomplished 
by interchanging the order of integration 'with respect to r 
and Td, which of course alters the limits of integration, 
giving 


<Pb- 


.=ssa tm -Jica. 

it V r a —r 2 


1 (VO di? , (VO dr 

v J*» dr D r ° Jr D -J(r u —r) {r—r D ) J> * ” 


or 


<?K~ 


2H M t .nV— +" 

V v r a —r 2 7 t dr D \ r a 

(61c) 


Equation (61c) may be integrated by parts to give the alter- 
native form 

■SE5 f*W T B 'ip, 

v tc (r u —r D )-yjri—r D v 

Thus, either by equation (61c) or equation (61d), <P B is 
evaluated by a line integral across the vortex field along the 
constant s a Mach line. 

The x and y components of the perturbation velocity may 
be obtained by differentiation of equation (61c). 
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If equation (73) is rewritten, with the use of the identity 

\dff_ d 2 If 


/_b_ 

\dr a ds„/ dr* 

there results 
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Integration by parts of the fourth term in the second member 
of equation (73a) leaves 


2|6 d^ 

*»(£-•) 

L 

21/ da: 

T-J(r„—r 2 ) (r 2 — s„) 
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(73b) 


(The series of equation (64) may be obtained from equations 
(73b) and (62) by successive integration of equation (73b) 
by parts.) In a similar manner, 

2 d -p s _ 2i/(0) Vr a — r a 

M dy jry (r m —r 2 ) (r a — s„) ~ir(r w —s a ) Vr 2 — 

2 pwafflY ■ » _ a a ■ JS1*> 
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( 74 ) 
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or 


dp H 

y(>+fe) 

m o) , 
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br D arD 
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(74a) 

Thus the perturbation-velocity components (associated with 
H) on the surface of the wing are given by simple line integrals 
across the vortex field along the constant s„ Mach line. For 
a given wing plan, the line integral of equation (73b) (and 
the first line integral of equation (74a)) is a function only of 
the shed vorticity and of the coordinate s a (as was similarly 
noted for equation (64)). The form of equation (73b) 
suggests that the equation might be more general than the 
example for which it was derived. 

If the Kutta-Joukowski condition is to be satisfied along 

djff 

the subsonic trailing edge, the value of that appears in 

equations (60a), (61a), (73), and (74) must satisfy the integral 
equation 


H(0) j p(»-) dr a dr ° U C 
Vr a — s w -y/r 2 —r D M Jt 


-q O' y) (Us 

bd 2 Vs»— S 


(75) 


where a depends on r 2 (s w ) and s. This expression results 
from adding the terms of equations (44) and (73b) that give 

infinite values of along the subsonic trailing edge (that is, 

those terms that contain the factor — 1 ^ V r »~ r *) 

and equating these terms to zero. An alternative evaluation 
of H may be obtained from equation (22a) and the known 

values of o or ^ for solutions satisfying the Kutta- 

Joukowski condition 

J'f d^e) * (76) 


If in equation (76), <P T and <P B are divided into those com- 
ponents associated with and independent of H (as in equation 
(58)), only the part associated with II remains after the 
subtraction. The part associated with II is zero from — <» 
up to the s u = 0 Mach line and has the value (from equations 
(73b) and! (75)) 
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(77) 


from the s w =0 Mach line to the subsonic trailing edge. 
(The + or the — sign is used for the top or bottom wing 
surface, respectively.) Application of the transformation 
equations (16) for. constant values of y (which infers that 
dr u =ds w ) to equations (76) and (77) then gives 


H =mf 0 


If (*»+fte 


ds u 


(del Q ■"» r*» 

^s a —My^r t J*i( r 
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(78) 


where x 2 represents the plan-boundary equation x=x 2 (y). 
This value of II may then be inserted in equations (60) and 
(74) to complete the evaluation of perturbation-velocity 
components associated with the discontinuity of sidewash 
behind subsonic (or supersonic) trailing edges. 

Acceleration-potential derivation of solutions satisfying 
Kutta-Joukowski condition. — Additional insight into the 
origin of. pressure forces may be gained by deriving the 
expressions for the influence of the upwash field in terms 


dp 


This deriva- 


of the acceleration potential ^ rather than <p. 

tion may proceed either from equation (42a) and figure 4 (a) 
or from equation (44) and figure 5 (b) under the assumption 
that X is finite along the wing boundary in conformance 
with the Kutta-Joukowski condition. Both derivations 
yield the same result, but the second approach is simpler 
to justify. If equation (44) is applied to the virtual wing 
region S Dt x of figure 5 (b) for both top and bottom surfaces, 
there results 
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(79) 


The last two terms arise in the evaluation of the surface integral from the discontinuity of streamline slopes along 
the trailing edge. Inasmuch as equation (79) holds for all values of r„, there results 
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(If equation (27a) is differentiated partially with respect to r and So, addition of the two resulting equations gives 
equation (79a). Inasmuch as equation (27a) excludes the effects of shed vorticity whereas equation (79a) does not, 



87 


AERODYNAMICS OP THIN FINITE WINGS AT SUPERSONIC SPEEDS 


the inference may be drawn that equation (27a) could include the effects of shed vorticity for time-independent 
cases by the addition of a function of r—s D . This fact was shown in another manner in reference 15.) With reference 

to figure 5 (c), ^ on the surface of the wing is 
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Elimination of A from equations (79a) and (SO) gives 
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Equation (80a) is the same as equation (65) except that the 
term generating infinite loads along the wing-plan boundary 
is missing. Thus, when the brace of equation (65) is zero, 
the Ivutta-Joukowski condition applies and equation (80a) 
results. The vorticity in the wake of the wing under such 
circumstances is of sufficient strength to cancel the contribu- 
tion to of the line integral along bd. 

Three approaches have been presented to show that the 
Kutta-Joukowski condition can be formally satisfied along 
subsonic trailing edges of thin wings at supersonic speeds. 
Each of these approaches implies a discontinuity in sidewash 
(that is, a vortex sheet) behind the wing of sufficient strength 
to cancel the term in equation (44) evaluated by the line 
'integral along bd. This term appears to depend only on 
those wing slopes lying on bd, which can only influence the 
flow field along or downstream of this Mach line. On the 
other hand, the contribution of the function 77, which opposes 
this integral for Kutta condition solutions (equation (75)), is 
determined by the discontinuity in sidewash established 
upstream of this Mach line (except for the isolated point d). 
The implication is that the Kutta-Joukowski condition 
cannot be satisfied without an interchange of cause and 
effect. Once the flow is established, however, the boundary 
conditions may be satisfied by solutions fulfilling the Ivutta- 
Joukowski condition. Application of the Kutta-Joukowski 
condition would seem to imply at least one of the following: 
(a) Downstream disturbances might be felt upstream through 
the wake of the wing or through the boundary layer whether 
initially or continuously. The apparent interchange of 
cause and effect in the linearized calculation would then be 
negated; (b) only the vortex line from the point d contributes 

9X8848—51 7 


to along the line s w (This argument seems to be in opposi- 
tion to equation (73b).); (c) the wing slopes must change 
continuously in approaching the line bd (as they would if 
they were represented by a power series) so that the apparent 
independence of the line integral of forward wing slopes 
might not be real. The application of the Kutta-Joukowski 
condition at supersonic speeds is still an open question, 
which should be determined by experiment. 

In order to illustrate the Kutta-Joukowski solutions, the 
perturbation-velocity components of a thin, flat, trapezoidal 
wing with a subsonic trailing edge on the tip may be derived. 
The wing equations are s— — hr and s=k 3 r. Either by inte- 
gration of equation (67) or when the singular term of equa- 
tion (48) is dropped, ^ may be derived as 


Uct (£i+l) __ t Ihihr^—Su) 

bx~0* V JfcW*. )s w 


(81) 


(Equations (78) and (73b) also lead to equation (81) for this 
example.) From equation (70), 
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This value may be inserted in equations (62) and (63). 
Addition of calculated from equation (63) or (74a) and 
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from equation (49) gives 
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Likewise, from equations (72) and (27e), 
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The values of the perturbation velocities, equations (81) to 
(83), are constant along radial lines from the origin showing 
that the flow is conical. This fact resulted from the use of 
straight-line plan boundaries, but was not otherwise assumed 
in the derivation of the perturbation-velocity components. 

Circulation and lift. — -As in subsonic-wing theory, the lift 
per unit span is directly related to the circulation, as shown 

in reference 25. The lift increment per unit span ^ is 

given as 

— 'ife-T?)* (84) 

where the subscript t indicates along the trailing edge. In- 
tegration of equation (84) gives 

3 £- -pU(<p B - PT ) t =2 P UH(y) = P UT (84a) 


where r is the circulation evaluated along a wing chord. 
Equation (84a) is the familiar subsonic .relation and may be 
used to evaluate experimental wing lifts from wake measure- 
ments, or to evaluate the theoretical lift of a wing as a line 
integral across the span. The velocity potential along the 
trailing edge is, of course, a function only of y, x having been 
eliminated by means of the trailing-edge boundary equation. 

Computed examples of perturbation-velocity compo- 
nents. — The theoretical expressions for the velocity potential 
near thin wings at supersonic speeds describe the flow char- 
acteristics. The nature of these flows may be clarified by 
numerical examples for flat-plate wings. The perturbation- 
velocity components in regions influenced by the tip 
of a trapezoidal wing are computed to illustrate the effects 
of subsonic leading and trailing edges. In addition, load 
distributions associated with angle of attack, uniform pitch, 
and uniform roll are illustrated for a complete wing. 

The effects of subsonic leading edges on the perturbation- 
velocity components may be illustrated for the trapezoidal 
wing of figure 13 (a). The supersonic and subsonic leading 

T 

edges are defined by the equations s, =— ^ and s s =2r. If 


these values are substituted into equations (48), (49), and 
(27e), the perturbation-velocity components in the 2=0 
plane are obtained. These components are presented in 
figure 14. 

The two addends that give ^ and are also included 

on figures 14 (a) and 14 (b). The arc-tangent expression 
controls the value of these two velocity components near the 
Mach line from the subsonic and supersonic leading-edge 
intersection, but the inverse square-root singularity becomes 
increasingly important near the wing tip. 

The value of ^ rapidly changes near the innermost Mach 
line from the origin, reaches an extreme, and then becomes 


unbounded near the tip. 




The quantity on the other 


hand, changes monotonically from a positive value on the 
innermost Mach fine to an unbounded negative value along 

the tip. The value of ^ (fig. 14 (c)) is likewise unbounded 

in the upwash field in the vicinity of the wing tip, but rapidly 

Qy 

drops as ~ increases. The upwash is zero along the outer- 
most Mach line from the origin. 


I 

r 




(a) Subsonic leading edge. 

(b) Subsonic trailing edge. 


Fioure 13.— Trapezoidal-wing plan forms. 
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(b) v component. 

Figure 14. — Perturbation velocities near tip of trapezoidai wins with subsonic leading* edge. 
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(c> z component. 


Figure 14.— Concluded, Perturbation velocities near tip of trapezoidal wing with subsonic leading edge. 


If the trapezoidal wing is flown backwards, the equations 
for the wing boundaries are «i=— 2r and s 2 =0.5r. In this 
case the wing tip is a subsonic trailing edge for which the 
Kutta-Joukowski condition, may be (arbitrarily) imposed. 
The three perturbation-velocity components, which were 
calculated by equations (81) to (83), are presented in fig- 
ure 15. — 

Both and ^ monotonically change from the two- 


da; 


dimensional regions to the wing tip. In the upwash field 
— -'<—^<0, ^ is zero and is constant, corresponding to 
a constant strength of vorticity in the wake of the wing. 
The velocity component is constant across the wing, then 
monotonically increases in the upwash field from the value 
along the wing tip to infinity along the line ~=0 (corre- 

1C 

spending to the limiting vortex line). 


Many of the equations presented heretofore are given in 
terms of the wing slopes cr. These slopes, however, are 
related to the induced z component of the perturbation 
velocity by equation (19). The formulations are cq-ually 
valid when w is induced by either steady wing motions or 
by wing geometry. The effective wing slopes on a flat-plate 
wing associated with lift, uniform roll, and uniform pitch 
are, for example, 

lift: 

<7 j> 

roll: 

( T jr 

pitch: 

(JTjf 

where m and n are the rates of roll and pitch and rjo and { 0 
are the coordinate distances to the roll and pitch axes, 
respectively. If these wing slopes are substituted in equa- 
tions such as equation (44), the load distributions associated 
with lift, steady roll, and steady pitch may bo determined. 
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( 85 ) 
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Pv 

X 

(c) z component. 

Figure 15.— Concluded. Perturbation velocities near tip of trapezoidal wing with subsonic trailing edge. 


The following table lists the type of equation that could 
be used to evaluate the velocity potential and the x compo- 
nent of the perturbation velocity for each of the wing 
regions of figure 16: 


"Wing 

region 

Equation for velocity 
potential 

Equation lor || (Kutla-Jouiowski condition 
imposed along subsonic trailing edges) 

I 

(17) - 


n 

(29) 

(44). _ 

hi 

(53) or (29) and (61b) 

(67) . . 

IV 

(3D 

Not presented. Equation is similar to (44). 

V 

(34) and (611)) 

Same as region IV, except term oi form 
(73b) from subsonic trailing edge is de- 
leted . 

VI 

(34) and (61b) applied to both 
subsonic trailing edges 

Same as region IV, except both terms of 
form (73b) from tiro snbsonie trailing 
odEHS are deleted ..... _ 


The calculations of reference 17 directly applied the equa- 
tions for to obtain the load distributions. 

The load distributions associated with lift, roll, and pitch 
for the wing of figure 16 are presented in figure 17. These 
distributions (reference 17) are shown in contour-map fashion 
giving lines of constant loading. The dashed lines are 
either the rolling or pitching uses or the Mach lines origi- 
nating at discontinuities on the wing-plan boundary. 

The wing loads associated with angle of attack (fig. 17 (a)) 
approach infinity along subsonic leading edges and are 
zero along subsonic trailing edges. The lifting pressure 
coefficients are conical in the forward central portion of the 
wing, generally decreasing toward the rear. The loadings 
are negative on the rearmost portions of the wing. This 
reversal of lift results from the pronounced effect of the 
upwash on the downstream parts of the wing. 
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As in the case of the lifting wing, the loads associated with 
uniform roll and pitch (figs. 17 (b) and 17 (c), respectively) 
are infinite along subsonic leading edges and zero along 
subsonic trailing edges (to satisfy the Kutta condition). 
The lines of constant load distributions associated with roll 
(fig. 17 (b)) are nearly parallel to the roll axis on the forward 
central portions of the wing. Considerable distortions are 
evident toward the rear, with pressure islands occurring. 
Positive pressure coefficients appear on the front and back 
parts of the wing, with negative values between. 



Figure 16.— Form of wing analyzed. 11= t /2. 


The load lines associated with pitch (fig. 17 (c)) are roughly 
parallel to the pitch axis for the central regions ahead of that 
axis. Pressure islands occur at the rear of the wing. An 
interesting observation is that the loadings associated with 
pitch change sign ahead of the pitch axis. 

The force distributions of figure 17 may represent the 
loadings of a variety of wing-plan boundaries. The wing 
plan may be altered or modified along any shaped supersonic 
trading edge without altering the loadings in any region for- 
ward of that edge. Likewise, if other axes of pitch or roll 
are desired, the new loadings may be obtained by super- 
position of values shown. If axes Vi and fi are desired, 
equations (85) become 
roll: 

ffr=— Jj in— Vi) = ~Jj (v—Vo) —Jj (Vo—Vi) 


pitch: 

—jj (i — &) = —jj (£— (?o — ?i) 

The shift in the axis (Va~Vi) thus superposes a lifting load 
distribution for a wing at angle of attack equivalent to 

(^o— *?i) on. the original loading associated with roll. 

m— TIME-DEPENDENT FLOWS 

The concepts that yield steady-state solutions for the 
velocity potential of finite wings at supersonic speeds may 
also be applied, in part, to time-dependent problems. The 
evaluation of the upwash field at a local point, however, is 
complicated by the fact that a knowledge of the wing-slope 
time histoiy is required. If the wing slopes do not vary too 
rapidly with time, approximate solutions may be obtained 
that nearly satisfy the boundary conditions in the plane of 
the wing. The theory is illustrated for a linear variation with 
time of the wing slopes, which gives an exact solution. 

Evaluation of isolated regions of upwash off subsonic 
leading edges. — Equation (23) may be written for the wing 
of figure 4 (a) as 

C r D Xg, 6 dr ds C t d P»z(0 (<t b — a- T ) dr ds . . 

Jo J*«( r ) - \l(r D —r ) (s D —s) Jo J»m 2 Y(r D — r) (s D —s) 

For steady-state solutions the slopes <r and X are inde- 
pendent of r D ; the integration of equation (86) with respect 
to s is therefore independent of r D and the reduction to the 
first-order integral equation (27a) is permitted. For the 
unsteady problem, the slopes <r and X of equation (86) are 
functions of the time delays r* and t„, which depend upon 
r D and s D , and the value of the integration with respect to s 
depends explicitly on the value of r D . The reduction of 
the second-order integral equation to a first-order equation 
by equating the integrations with respect to s is therefore 
seldom justified. No published closed-form solutions of the 
second-order integral equation (equation (86)) are known 
for arbitrary wing slopes, but approximate solutions are 
obtainable. 

If the wing slopes generally are evenly distributed, those 
slopes near the forward Mach lines (represented by 
(r D —r) (sjj — s)=0) will carry the most weight in the evalua- 
tion of the second member of equation (86) . Near these Mach 
lines, however, the nonlinear time delays of equation (18) 
become nearly linear with respect to either r or s and the 
two time delays r„ and r t are nearly equal. This fact 
suggests that the slopes appearing in equations (86) may 
be expanded in a power series in terms of the variable 
2 

f~M(3c V( r D — r ) (sd~s). The subscript a' may refer to the 

slopes or the derivatives of the slopes evaluated at /= 0, 
which is the same as applying a time delay of 


r D — r , s 
r “' j8c + j9e 


(87) 
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(The notation r a > will also be employed for point (r a ,s„) 
rather than ( r D ,s D ).) The expansions for a and X for the 
two time delays r a and n have opposite signs attached to the 
odd power terms in/. The odd terms therefore disappear 
when the two series are added 

(eg — o't), i.s =2 ^(ffs — o>)a'+ 2 | — ~ B ^)f- T ^ a 

css) 


X..,=2(x..+I • ■) 


(89) 


If the wing slopes do not vary too rapidly with time, the 
higher-order terms of equations (88) and (89) are negligible 
compared with (<rg— a T )a r and X tt '. This approximation is 
equivalent to assuming that the slopes change linearly with 
time in the interval r a — r t =2/. This approximation simpli- 
fies equation (86) to 


C t d T j d X g >dr ds C r D C 

Jo r) (s D — s) ~ Jo J* 


‘•sW (a B — o>) a ’dr ds 
*i(o 2 V {r D —r)(s D — s) 


(90) 


A solution of equation (90) that in form resembles equation 
(27d) is 

Mr,s,t)= n } dp ( 91) 

2 xyS— 8 S ($— V) 


where the wing slopes (trg— <r T ) of equation (91) are evaluated 
at times t—r a " and 

s—d 


/3c 


(92) 


In the analysis that follows, the approximate value of X 
from equation (91) (which is the exact solution of equation 
(90)) is applied as if it satisfied equation (86). 

Aerodynamic effects of isolated subsonic leading edge. — 
Analogous to the steady-state equation (28), the contribu- 
tion <p D of the upwash field S B to the velocity potential on 
the surface of the wing in figure 4 (b) may be computed by 
equations (17) and (19) as 


<?D= 


U dr Xg.a ds D . . 

2A/tJo /(r„— r) J»jcd 4^—Sd 


The streamline slopes X afr are evaluated for the point 
(r, Sd) at times ( t—r a ) and (f— r b ). On the other hand, the 
evaluation of X in terms of wing slopes by the approximate 

g g 

equation (91) includes an additional time delay of 
The combined time delays become 


/Sc 


®I£ ^ I r H T t ^ l-f T-. T 

/Sc /3c r ) ( s « — s a) 

Sm — S , r w — r 2 ^ r i 

/3c ^ /Sc r ^ Sa J 


(94) 


Substitution of equation (91) in equation (93) then yields 


U (* r j(*») dr ds D C'iW Vs2 — s(tr B — a- T )~_ids 

rD 4k 3 !? Jo ^{r w —r) J's( r ) V(s D — s 2 ) (s a —s B ) J»i( r ) s B —s 


(95) 


The order of integration with respect to s and s D may be interchanged in equation (95), so that the potential <p D may be 
represented as an area integration on the surface of the wing: 


U 

4 


fM*-) dr /V») , , 

Jo V(r B — : r) J«iW S2 8 J*a( r ) (s D - 


(cr g O’ t) 7. i dSp 


-s) / (Sp Sj) (s w S b ) 


(95a) 


If the approximate equation (95a) is combined with the con- 
tribution to the velocity potential on the top wing surface 
(exclusive of the upwash field), as calculated by equations 
(17) and (19), there results 


= u r r M^drds u_ r r 

V ~ 2l75r J Js„,i V {Ta—r){s a —s) 2 Mt J Js 


Q dr ds 


where 
Q= (ff r ) a , s 4 


V (s a — s) (s 2 —s) 


r*» (o-g— g>); 

J* J (s D — s) V (s D — 


s »,j V (r w —r)(s w —s) 
(96) 


, ds B 


w $Z>) 

(97) 


(Because X becomes unbounded along the subsonic lead- 
ing edge, the exact time delay is employed in equation (93) 
even though X is given only approximately by equation (91). 
This choice also matches equation (96) to the exact solution 
of Garrick and Rubinow (reference 11) along the s a =0 Mach 
line. For many problems involving slowly changing wing 
configurations, (o- r ) a ,» and (<rg— cy)- r of equations (96) and 


(97) may be expanded in series form similar to equations 
(88) and (89) 

(c'r)a,s == 2((T2’) a '-r • • ■ 

(eg — or)- y= 2 (erg — <fr)a'~ f~ • • • 

The value of Q then becomes Q=(o-g+c- r ) a . and the velocity 
potential, equation (96), simplifies to 




U 


XX 


(<T T )a' dr ds U 


Af 7T J Js„,i -y/(f w — r)(s m — s) Mt 


/X 


(°~a~{~g~r)c' dr ds 
s -,i 2 V (r„— r) ( s„— s ) 
(96a) 


Equation (96a) illustrates a theorem that for moderate 
time-dependent wing motions, steady-state equations (see 
equation (29)) may be employed if the steady-state wing 
slopes are replaced by time-dependent wing slopes utilizing 
a time delay of 

_ r a — r , s a —s 
r *'~ /3c + /3c 

Equation (96a) gives exact solutions (linearized theory) if 
the wing slopes vary linearly with time.) 
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A simple example in which the angle of attack of all wing elements of a flat plate varies linearly with time (reference 18) 
illustrates the equation. The wing slopes may be expressed as 


ctb = a "F rut— — ctt 


(98) 


These effective wing slopes may correspond to a constant acceleration of mU in the z direction of a wing at angle of 
attack a. Substitution of equation (98) in equations (96) and (97) (or equation (96a)) and conducting the manipulations 
for a wing whose leading-edge equation is Si= — hr give 


v “m\[ a+mt ~Wc { ' k '~ l)r ‘~Wc (f, - !) ( 1_ l) r »] V<r.-rO(*i+iiri)+ 

Substitution of equation (99) in equation (2a) evaluates the pressure coefficient as 


(98) 




( 100 ) 


where r 2 is evaluated at s=s„. The steady-state solution (equation (45)) results if m=0. The load distributions of a 
family of wing-plan boundaries is obtained by choice of the equation r=r 2 (s„). The solution for the infinite swept wing 
results when s to =r 2 (s IO )=0. Inasmuch as the higher-order terms of equation (88) are zero for this example, the solutions 
given by equations (99) and (100) are exact. 

Calculations of perturbation-velocity components. — As in the time-independent cases, expressions may be derived from 
the approximate equation (96) for the perturbation-velocity components. These components are obtained in a manner 
similar to the derivations of equations (44) and (45), which are presented in appendix B. Direct computations of the 
perturbation-velocity components in oblique coordinates yield 


dta 


dr w 


&} I 


57 ** 


u 


I s * V(n»— r) (s„— s) 

U C [Qrr )o-q ds 
hd V (r«,— r) (s „— & 


2 Mr. 


d(°>)<..i> t j 

n* _ dr ar as 
s ‘ Vfo,— r) (s» s) 


U 


X 


fffrla.s dr? 


u 


22 / 71 - Jab ^(rw—r) (s w —s) 22 /xJb 0 d V(r„— r) (s„— s) 


X 


Q ds 




U 


XX 


-s) 


^ dr ds 
ds 


( 101 ) 


U 


ds„ 22 /tJ Js, V (r w — r) (s„— s) 2Mt, 


a (°~r)«.i j j 

r r ds ar ds 

J Js V(t»— r) (s„— s 


E7 f —U T )a.t,dr 


ILC 


-Q dr 


s) 22/7rJab V( r »— r) (s„— s) 22/vJbod V(r K — r) (s„— s) 


22/71- J b 


bd r ) (S«~S) 


( 102 ) 


where the differentiations with respect to r and s are conducted with (r ra — r) and (s„— s) held constant. 

The evaluation of Q is especially simple for the line integrals along bd, because only the point d is common to the line 
bd and the upwash field so that s„«s D for the evaluation of the wing slopes. The quantity (<r B — v r )a,&is then independent of 
s D and can be taken outside of the integral in equation (97) to give 


Q — (c t) a.,b~T (vs — Vj-ls' (for So ~S a ) 

Syr - "S , r %a r 2 (s,ff) 

TS '~~pT' t pc 


(103) 


where 


( 104 ) 
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Equations (101) and (102) may then be written 

S (2 7 7 b(<Ty)c»fe J J 

u c C 5r dr “ s u rr ar dr ds u_ r ds v_r Qds 

£>r„ 2 MrJ Js, V0\s—r) (s„— s) 2iLfrJ Jsj y(r„— r) (s K —s) Jib V(r„— r)(s„— s) 2iU*-Jbod V (r u — r) (s„— s) 

£ f — (ioia) 

2'Mir\r w — r 2 (s„) Jbd v's»— 5 


and 


— dr ds d(or 

a*_ jj_ r r as u_ c r a 

as* 2d/ir J Js t V (r„ — r)(s„— s) 2ALr J J S) V(r„— r) 



u 

C — (o'rlg.jdr 

U 

C ~Qdr 

U JT a 

C ( or B —cr T )a'ds 

2Mt , 

V(r„— r) (s»— s) 

2A/t. 

Jsci V(r„— r) (s„— s) 

2Miry'r u —r 2 , 

Jbi yS w —S 


The a? component of the perturbation velocity is calculated from equation (15) 

a^ Af / a ip . a y \ 

ar"2/3 Var,,, 5s 


(102a) 


(105) 


When equations (2a), (105), (101a), (102a) , and (96) are combined, there results for the pressure coefficient on the top wing 
surface 




(*Q 

V ar' 


d<3, 2/3 

' a.s ^ mu 


drds 


V( r w —r ) (s w — s) 


2/3y 


J/J 


( 


d(°~r)g.& , a(o~r)g.s t 2ft a(o~r)g.i> ~| j j 

^ "T" -V IS f TT ^ r I (*/ UO 


ar 


a s 




^(r a —r) (s u —s) 


1 _ r (o'rlg.afds — dr) 1_ j' Q(ds-dr) 

2/3 TT J fl4 y (jr a — r) (s„— s) + 2/35r Jw<f V (r„ — r)(s w — s) 



p 


(q- b — er r )t' ds 
Vsg>-s 


(106) 


where the partials with respect to r and s are evaluated with (r„— r) and (s„— s) held constant. One way to perform this 
differentiation without error is to replace ^ with (r u — r), s, (s„ — s), and t held constant by with all of the 

variables except the one of differentiation held constant, and so forth. With this innovation, the operator 


of equation (104) becomes (^r+W+^+^+MJ bt)‘ 

The approximate solution for the pressure coefficient that satisfies the Kutta-Jonkowski condition along subsonic 
trailing edges is equation (106), except that the integral along the line bd is deleted. This conjecture has been 
demonstrated for a linear time variation of wing slopes in reference 26. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, June 17, 1949. 



APPENDIX A 


The following symbols are used throughout this report : 
a 0 , a u a 2 , . . . a n coefficients of power series expansion of 
H(y)=a a J r a l y+a 2 y i . . . 


SYMBOLS 


c 

Up 

C 

F 

F n 

f 

/ 

9 

9(s„) 

H 

k 

L 

l 

M 

m 

n 

V 

Q 


R 
R' 
r, s 

r, s 

S 

t 

t' 

u 

v, v 
w 


A'(/v, s) 

1 ( r , Sp.) 

V, z 

x', y', z' 


pressure coefficient, ?- — — 
t 2 p*U* 

speed of sound 

suction force in flight direction 
suction force normal to wing edge 
function 

Mfc V(ri>—r) 

integration function of time (herein con- 
sidered as constant! 

Ua 

times bracketed portion of equation (66) 

integration function of x-Ut and y 

constant- 

lift 

distance along plan boundary 
free-stream Mach number 
rate of roll or rate of change in angle of 
attack 

rate of pitch, summation index, or coordi- 
nate normal to curve 

static pressure or function -y'rj— s w 

(o' r) <J, l 

>i(n (ctb — o'r);.T dsp 


V(s E -s) (s a -sK ,, 


p 

J»1 


(r) (Su— s) V (Sd—Sz) (s„— Sd) 

= V(z-g) 2 -ff 2 G /-i?) 2 -iS 2 (z-r) 2 

= ■\ix' 2 +y' 2 +z' 2 

oblique coordinates whose axes lie parallel 
to Mach lines in 2=0 plane 
point of tangency of wing plan boundary 
and foremost Mach line 
plan-form area 
time 

= u-mt+^Y 

free-stream velocity 
variables of integration 
2 component of perturbation velocity meas- 
ured positively outward from surface of 


Cartesian coordinates (free stream parallel 
to r-axis) (The x- and y - axes coincide for 
top and bottom wing surface. The z-axis 
is measured positively outward from 
either surface of wing.) 
transformed Cartesian coordinates 

(x f =x, y' = Py, z' = VT^H 2 z) 


a 

0 

r 

x 


p 

a 


T 

<P 


Subscripts: 

0 

1 , 2 , . . . 

a, b, a', a", a, 
b, and a' 

B 

D 

H 

T 

t 

w 

Examples: 

r 2 

Si 

Sua+i) 


X(^& 

X 0 

X^ 

(ffs — crr)a,b 


<Pd 


T(re+1) = 


Cartesian coordinates in x, y, z directions, 
respectively 
angle of attack 

cotangent of Mach angle, 1 

gamma function or circulation 
slopes of stre amlin es (measured in incon- 
stant planes) in z=0 plane between wing 

V) 

boundary and foremost Mach line, X=^j 
density 

effective wing-section slopes measured in 
i?= constant planes, <r=^ 
time delay 

perturbation-velocity potential (variation 
from free-stream velocity) 

free stream, axes of roll or pitch, or without 
vorticity 

numbered areas or wing-plan-boundary 
equations 

time delays t,, rj, r,, tj, and t,', 

respectively 
bottom of wing 
upwash field 

associated with shed vorticity in upwash 
field 

top of wing 
trailing edge 
'ring 

curve r=r 2 (s) 
curve s=s t (r) 
wing area 1 plus 2 

derivative of curve r=r 3 (s„) with respect 
to s„ 

slope at time t—r a plus slope at time f— r# 
part of X remaining if H is zero 
part of X associated with H 
difference between bottom and top wing 
slopes at time t— r„ plus this difference 
at time t— t„ 

contribution of upwash field to velocity 
potential on top wing surface 
n\ 

M^Su—s+Ty—r) +2 ^(r w — r) (s a —s 


wing 


Mfic 

1 C r w w dr 


M(s u —s+r u —r)—2\(r u —r) ( s u —s ) 

Air \r m —r 


- M0c 

1 f*« w ds 


r D —r , s D —s 

J—cb \S w S 

Ttt ' 

1 

1? 


Ta" 


S D — S 

J3c 

Hf(g a — s+r M — r)+2^(r u — r)(s„— s D ) 

M(3c 

M(j} w —8+r u —r) — 2-J(,r u —r) (s„— s p ) 


Su S . Ty; T 2 [S jp) 


Mj3c 


0C 


0C 
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APPENDIX B 

DERIVATION OF EQUATIONS (44) AND (45) 


The derivation of equation (44) in Cartesian coordinates is 
presented in reference 15. An alternative derivation in 
oblique coordinates is presented herein. The r„ and s w 
components of the perturbation velocities are obtained for 
each of the areas S u ,i and S ait of figure 4 (b) in a manner 
analogous to the derivation of equation (42) of the text. 
The velocity potential at point (r a , s a ) associated with 
the area S w ,i is 


b 



U f C cr T dr ds _ 
(Pl Mr J Js„,i -y/(r w —r) (s w s) 


(Bl) 


The potential an infini tesimal distance in the r direction from 
the point (r„, s„) may be written 

.a*,,. t/ r r (*'+w ir -) drds 

r. Mr] Js Wil yl(r a -r) (s u -s) 

Udr„ C <r T ds Udr a C 

Mt Jab V(r w — r) (s m — s) Air Jb 


m ^(r a —r) (s w —s) 


Subtraction of equation (Bl) from equation (B2) leaves 

dr (Jo 

a g i = u r r dr aras 

dr„ A/xJ -J (r w —r) (s w —s) 

TJ or ds TJ p or ds 

Air Jab -}(r w — r) (s w s) Aixjbd -J (r u —r) (s„— s) { 


In a similar manner 


(J r (J~ 

U CC ds araS . TJ r <r T dr 
ds a StJ ^(r*— r) (s B — s) + AfxJ a b V(r B — : r) (s a —. 


f 


ffr( 


U_ dri 

Mt ds v J bd V(r B — r) (s w —s) 


s) 

(B4) 


The velocity potential associated with the area S a _i (and 
the upwash field) follows from equation (29) as 


TJ C C C o b 4“ o J) dr ds 

^ AiirJ 2 V (r„— r) ( s a —s ) 


m 


b 



The potential at (r„+dr«,,Sa,) may be written as 


<p 2 + lrl dr ”- SJk 


s i i d(og+0r) 
cr B +0-7.4 ^ ar v 


drds 


2 V(r„— r)(s„— s) 
( <7 B 4 & t) dS 


Udr w r 

Mir Jbod 2 y ( r w —r ) (s„— s) 


Udr a f 
Mt Jbc 


( 05 + 0 - 7 .) ds 


CB6) 


/bd 2 V(r»— r) ( s„—s ) 

Subtraction of equation (B5) from equation (B6) leaves 
1 d(g- 8 +<r r ) . 7 

u_ r r 2 __dr aras 

dr„ Mt} Js»,j V(r„ 


r) (s a —s) 

TJ C ( &B + d~r ) ds ■ U r (og~h & r)ds 

A/xJbod 2ii(r a —r) (s„—s) ‘ Aixjbd 2 V (r„— r) (s„—s 


Similarly, 

1 d(os+o>) , , 

d^ = u_ r r 2 ds drds 

ds„ Mt} J Sip ,2 V (**»— r) (s»— s) 

(05 + 07 .) dr f/ dr a 


(B 7 ) 


JLC „ 

Aix J bod 2 V (r„— r) (s„— s) A/x ds, 


L C (og+o^ds 
»Jbd 2 V(r„— r) (s„— s) 

CBS) 


Addition of equations (B3) and (B7), and also equations 
(B4) and (BS), gives the perturbation-velocity components 
in the r and $ directions, respectively, as 

drds 

d® _ fj r r dr a 

dr„ A/xJ Js».! V(r B — r)(a„— s) 

1 d(o5+0r) ,7 

u r r 2 dr U_ r <r T ds 

Mt} Js»,j ^(r w —r) (s„—s) A/xJ tt b -J(r w — r) ($ w — s) 


TJ C (og+ 07 .)ds £/ f (or 05 ) ds 

A/x J bod 2 V (r«.— r) (s K — s) A/xJbd2V(r«— r) ($ u —s) 


U r (or- 


(B9) 
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_U_ 

Mir 


C )<7 r r 

bs 


drds 


Sw f \ ^ {t uf ig s) 


_U_ 

Mir 


2 os 


U 


s. 


— ( 7 x dv 


s*.i V(r«,— r) (s w — s) 

U C — (<xs or) dr | U dr 2 C 

MirJbOd 2sj{r w —r) (s a — s) Mir Jl^Jbd. 2 V (r w —r) (s„— s) 


Mir Jab V (Tiff — r) (s a —s) 
(a r —<r B )ds 


(BIO) 


The x and y components of the perturbation velocity, which 
are given by equations (44) and (45), may be obtained from 
equations (B9) and (BlO) and the transformation equations 


btp M / bip . bp \ 

bx~2& \5r„ dsa,/ 

0 


b<p M / 

by~ 2 \c)Sb, 


bp ' 

br„ 


(Bll) 
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